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^ ■ Abstract 
■ 

' The A^/-flavour Schwinger Model on a finite space < < L and subject to bag-type 

boundary-conditions at = und = L is solved at finite temperature T = 1/13. The 
boundary conditions depend on a real parameter 6 and break the axial flavour symmetry. 
We argue that this approach is more appropriate to study the broken phases than intro- 
ducing small quark masses, since all calculations can be performed analytically. 
In the imaginary time formalism we determine the thermal correlators for the fermion- 
fields and the determinant of the Dirac-operator in arbitrary background gauge-fields. We 
show that the boundary conditions induce a CP-odd 6'-term in the effective action. 
The chiral condensate, and in particular its T- and L- dependence, is calculated for Nf 
fermions. It is seen to depend on the order in which the two lengths (3 = 1/T and L are 
sent to infinity. 



1 Introduction 



Over the past decades the Schwinger model ^ has proved to be an excellent laboratory for 
field theory because it turned out to shed some light on a couple of questions which naturally 
arise in realistic gauge-field theories, but lead to immense difficulties as soon as one tries to 
attack them directly. Longstanding problems of this type are the wellknown C/(l)^-problem 



1 30 1, the question whether QCD in the chiral limit shows a spontaneous breakdown of the 
chiral symmetry and the question about the nature of the chiral phase transition at ~ 200 
MeV 1,0. 

The Schwinger model is known to be the most simple model field theory which exhibits 
chiral symmetry breaking. However the quantization on the plane suffers from the deficit, 
that a naive calculation of the condensates {ipip) and (V'TsV') gives zero results and the correct 
values can be derived only a posteriori by using the clustering theorem Q. Whenever a 
symmetry is expected to be broken it is most recommendable to break it explicitly and to 
try to determine how the system behaves in the limit when the external trigger is softly 
removed. Thus it is most natural for both the Schwinger model and QCD to break explicitly 
the axial flavour symmetry and to investigate how observables do behave in the limit where 
the symmetry is restored. 

The most direct way to do this is to introduce small fermion masses and to try to determine 
how the chiral correlators behave in the limit where these fermion masses tend to be negligible 
as compared to the intrinsic energy scale of the gauge-interaction. Once these calculations 
do predict nonvanishing chiral condensates in the thermodynamic limit one can be sure that 
a spontaneous breaking of the axial flavour symmetry SU{Nf)A really takes place - however 
this is not a conditio sine qua non. There is however a technical obstacle to this approach: the 
value of the chiral condensates is related to the mean level density of the eigenvalues of the 
Dirac operator 0] in the infrared. Unfortunately, the spectral density of the massive Dirac 
operator is only known for very special background gauge-fields. 

In this paper we shall break the chiral symmetry explicitly by boundary conditions for 
the fermions instead of giving them a small mass. Although this version seems at first sight 
less natural, it has many advantages - both conceptual and calculational in nature. The 
most important point is clearly the fact that it allows for an entirely analytical treatment. 
In a previous paper we investigated QCD-type theories with A'^^ massless flavours on an 
even-dimensional {d = 2n) euclidean manifold M with boundary dM on which the boundary 
conditions studied by Hrasko and Balog Q have been applied. These chirality-breaking- 
(CB-) boundary-conditions relate the different spin components of one flavour on dM and 
are neutral with respect to vector-flavour transformations - so that the (gauge-invariant) 
fermionic determinant is the same for all flavours. For a simply connected M, e.g. a ball, the 
instanton number, which in four dimensions takes the form 

-^iLu ^uu d X, (1) 



327r2 



is not quantized and may take any real value Contrary to the situation on a compact 
manifold without boundary, on which q is integer [^], the configuration space is topologically 
trivial (i.e. without disconnected instanton sectors) H]. In addition there are no fermionic 
zero modes [^, 11 1 which usually tend to complicate the quantization considerably Q. 



Our previously cited work focused on the euclidean Aj-flavour U (Nc)- or SU (Nc)- gauge- 
theories inside 2n-dimensional balls of radius R. We computed that part of the effective 



1 



action reflecting the interaction of the particles with the boundary S'^~^. Here we investigate 
whether the approach of breaking the SU{Nf)A symmetry by boundary conditions can be 
extended to gauge-systems in thermal equilibrium states. In the imaginary time formalism 
spacetime is then a cylindrical manifolds M = [0, /?] x {space} and (anti)periodic boundary 
conditions for the (fermi)bose-fields in the euclidean time with period /? = l/T are imposed. 

Note that at finite temperature it is only the boundary of space and not of space-time 
where chirality is broken and it is a priori an open question whether this is sufficient to trigger a 
chiral symmetry breaking even in the one fiavour case. In addition there is a technical obstacle 
to extending the CB-boundary-condition approach to non-simply connected manifolds, e.g a 
cylinder. On cylinders the standard decomposition for the Dirac operator, on which the 
analytic treatment heavily relics, must be modified. The present paper is a technical one - 
mostly devoted to show how this difficulty can be overcome. From the physical point of view 
it is our aim to investigate how the breakdown of the chiral symmetry - when triggered by 
boundary conditions - in the one- and the multi-flavour cases is affected by flnite temperature 
effects. 

Here we shall quantize the Schwinger Model with action 

S[A,'iPI^]=Sb[A]+Sf[A,4^1^] 

Nf (2) 
M n=lM 



on the manifold 



M = [0, /?] X [0, L] 9 ix^,x^) (3) 



with volume V = PL. At flnite temperature the flelds A and are periodic and antiperiodic 
in euclidean time with period /3 and hence = and x^ = (3 are identified. This means that 
[0,/?] X [0, L] is a cyhnder with circumference (3 = 1/T. At the spatial ends of the cylinder 
(i.e. at = and x^ = L) specific CB-boundary-conditions are applied. Then there are no 
fermionic zero modes (see next section) and the generating functional for the fermions in a 
given gauge-field background A is given by the textbook formula 

ZF[A,r]\v] = det{ip) e/''^^^^)"''?. (4) 

We shall see that these CB-boundary-conditions indeed generate chiral condensates for 
any finite length L of the cylinder. However in the limit /? ^ oo,L ^ oo the condensates 
will only survive for the one flavour case and this only if the limit /? — > oo is taken before the 
limit L ^ oo. 

During the calculations the following abbreviations are used for notational simplicity: 

This paper is organized as follows : In section 2 we discuss the CB-boundary-conditions 
to be applied together with some immediate consequences for the spectrum of the Dirac 
operator Section 3 is devoted to the question of how to decompose an arbitrary gauge- 
field on a cylinder. In section 4 we compute the fermionic Green's function with respect to 
CB-boundary-conditions in arbitrary external fields. In section 5 we determine the effective 
action after the fermions have been integrated out. Using the results of the two previous steps 
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the chiral condensates are calculated in section 6. In section 7 we show that the value of the 
chiral condensate crucially depends not only on the number of flavours but also on the order 
in which the two limits /? — > cxd and L — > cxd are performed. Finally we compare our result 
with the condensate generated by fractons on a torus of identical size and with analogous 
results of noncommutativity of the limits m — > 0, L ^ oo in the the usual small-quark-mass 
approach. In the appendices we derive the boundary Seeley-DeWitt coefficient used in the 
body of the paper. 



2 Chirality Breaking Boundary Conditions 

In this section we shall shortly review the boundary conditions as discussed by Hrasko and 
Balog together with their most important consequences H, 26]. 

Since Zp should be real we want iJp to be symmetric under the scalar product 



from which we get the condition 



M 



(X, ipi^) - {ilpx. ^)=i j xhni' = . (6) 

dM 

Imposing local linear boundary conditions which ensure this requirement amounts to have 
X^ln"^ = on dM for each pair, which is achieved by 

ijj = B%l) on dM with B^nB = -^in , 5^ = 1, (7) 

where 7„ = (7, n) = n^7^ = ^ and is the outward oriented normal vectorfield on dM. 
We shall choose the one-parametric family of boundary operators |Q] 

B = Be:= i75e'^^7n (8) 

which is understood to act as the identity in flavour space. These CB-boundary-conditions 
break the 75 invariance of the theory, making the Nf flavour theory invariant under SU {Nf)v 
instead of SU[Nf)L x SU{Nf)ji. Later they will be supplemented by suitable boundary 
conditions for the gauge-field. These boundary conditions imply that there is no net U{1)- 
current leaking through the boundary, since n ■ j = ip^^ni^ = on dM. 

In the following we shall make use of a Feynman Hellmann |16| boundary formula, which 
may be derived from (^,0J^) 

d % f 

-^>^k = 2 J ^kil ■ ^h^i^k = -Afc(V'fc, 75V'fc) , (9) 

where the denote the eigenvalues of ip. 

We choose the chiral representation 70 = o"i,7i = ct2 and 75 = o"3 . Then the boundary 
operators at the two ends of the cylinder read 

BL = -(^e Q)(atxi = 0) and BR = + (^g "J ] (at = L) . (10) 
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The most important properties of these boundary conditions are summarized as fohows ^j: 

(1) The Dirac operator has a purely discrete real spectrum which is not symmetric with respect 
to zero. 

(2) The Dirac operator has no zero modes. 

(3) The instanton number q = ^ J e^yF^y = J E is not quantized . The second property 
allows us to calculate expectation values of gauge-invariant operators as 

(O) = J{0)^ dfielA], where (11) 

dM^] = ^ dete(i^) e-^^[^l (12) 

Here D[A] is assumed to contain the gauge- fixing factor including the corresponding Fadeev- 
Popov determinant and {0)a denotes the expectation value of O in a fixed background 

Throughout 6 is the free parameter in boundary operators ([To|). We shall see that the 9- 
dependence of the fermionic determinant detg{i]p) can be calculated analytically. 



3 Decomposition and Deformation techniques 

In this section we present the decomposition and deformation techniques needed to determine 
the functional determinant of the Dirac operator on the cylinder with CB-boundary-conditions 
as given by (|7|) and (|10|) . 

On simply connected regions we have the decomposition eA^ = —e^ydu4> + d^x such that 
cFqi = A(j). On the cylinder there is a one to one correspondence between </> and eFoi if (j) 
obeys Dirichlet boundary conditions at the two ends of the cylinder. But cylinders are not 
simply connected, 7ri(M) = Z, and as a result the Polyakov-loop operators 

i.e. e Jq Aq mod 27r, are gauge- invariant. On the other hand, using the /3-periodicity of 
X and the Dirichlet boundary conditions on (p, the above decomposition would imply that 
Jq Jq Aq = 0, a condition which does not hold in general (take a constant Aq). This simple 
observation already indicates, that the correct decomposition of on the cylinder reads 

eAo= -di^ + dox 
eAi= +do(t> + dix, 

where (j) obeys Dirichlet boundary conditions at = 0, L and x fulfills x(0) + x{L) = and 
c E [0, 1[ is the constant harmonic part. To prove ( [l^ one Fourier decomposes the various 
fields and carefully handles the zero-modes of the Laplacian. The harmonic part can then be 
reconstructed from its values on the boundaries. 

The Dirac operator i]p= i^^{d^ — ieA^) may be factorized according to 

ilp=G'^i%G, where i^o = 7°(«5o + 27rc//3) + 7^i9i, (15) 
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and 

G=(^7 °)' 9:^e-^^^^^\ (16) 

The prepotential g is an element of the complexified gauge-group U{1)* = x Now we 
deform the prepotential and Dirac operator as 

.= e-'^i'P+ix) and ip^ = Clilpo Go, (17) 

such the deformed operator interpolates between the free and full ones: iJpa=i = "ip and 
iI/)a=o = ilpo- By using 

= -G^H , H=(~l^* M = -</)75 + ixl ■ (18) 



one finds for the a-variation of the integrated heatkernel of {ilpa 



\2 



-^(tr{e*-^}) = 2ttv{e''^-{H + H^){i]Paf} = 2t4( tr {e*^^"^'2</.75}) (19) 
da at 

and this formula will prove to be useful in section 5. 



4 Fermionic Propagator w.r.t. Boundary Conditions 



In order to calculate the condensates we need the Green's function Sg of the Dirac operator 
iJp on the cylinder subject to the CB-boundary conditions. This Green's function obeys 



{i]pSg){x,y) 
{BRSe){x\x^=L,y) 



6{x -y) 
-Se{x,y) 
Se{x^,x^ = 0,y) 
Se{x°,x^ = L,y) 



(20) 
(21) 
(22) 
(23) 



plus the adjoint relations with respect to y. The dependence of the gauge-potential has not 
been made explicit and the boundary operator -B^/zj is the one defined in (|l^). From the 
factorization property (|l5|) for the Dirac operator it follows at once, that is related to the 
Green's function 8$ of i]pQ as 



Se{x,y) = G-\x)Se{x,y)G^-Hy) 



(24) 



Indeed, since the field (j) obeys Dirichlet boundary conditions at the ends of the cylinder, g is 
unitary there and the boundary conditions (21-2^) transform into the identical ones for 



Se{x,y) 



S+- 



where the indices refer to chirality. 

The free Green's function on the cylinder of infinite length 
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Stherix,y) = — J2i-^T^ 











n/3 
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where 



-y^, is purely off-diagonal and thus chirality preserving, as expected, since ther- 



mal boundary conditions are chirality-neutral. To implement the chirality-breaking boundary 
conditions at the ends of the cylinder on can either augment Sther by pieces built from the zero 
modes (which themselves cannot obey the L/R conditions simultaneously) or by exploiting 
analyticity arguments. In either case we end up with 



I 

2^ 



E (-1 



Am+n) _ 27ric(g0//3-n) 



6 /fnm 
~ (1/ ^nm) 



(1/ ^nm) 
6 l°^nm 



(25) 



where r^m = ^ ^ir\ — {n(3 + 2imL) and Snm is the same expression with r] = x"^ + replaced 
by Prom this explicit expression one sees at once that the off-diagonal elements only 
depend on x^—y^ and become singular for x ^ y , whereas the diagonal elements depend on 
both x^ and y^ separately but are regular at coinciding points inside the cylinder. The sum 
over m respectively n in (25) can be performed by using 



E(-i 

z 

with the results 



ITT 



m + ia 



sinh OTT 



(— vr < X < vr) 



(26) 



Se{x,y) 



AL 



•E(-i 

z 



^—2ninc 



or 



sinh(7rr„o/2L) 
1 

sinh(7rs„o/2i) 



Se{x,y) = '-^^-Y.i-^) 

z 



e e 



sin(7rrom//3) 

^-2lTicSom/l3 



sinh(7rs„o/2L) 
sinh(7rf„o/2L) / 

-27ricsom//3 \ 



(27) 



sin(7rsom//3) 

Q-6^-2mcrom/P 

sin(7rsom//3) sin(7rfom//3) ^ 



(28) 



valid for c G [— ^, For calculating the chiral condensates we shall need the ++ and 
elements at coinciding points inside the cylinder. From (p^) we find the expression 



Se{x,x] 



±— y(- 

n£Z 



„±2innc 



sin(7r[^ — inr]) 



(29) 



which rapidly converges for low temperature, and from (|2^) the alternative form 



Se{x,x] 



±— y (-)" 

2/3 ^ 



mez 



g±27rc(^+m)/T 

sinh(7r[^ + m]/r) 



(30) 



which is adequate for high temperature. 

With (24) we end up with the following expressions for chirality violating entries of the 
fermionic Green's function on the diagonal 



5e(x;x)^ =eT2<^(")5e(x;x) 



(31) 

The free Green's functions S±± have been computed in (|29|) and (|30|). They depend only on 
the harmonic part c in the decomposition (14) of the gauge-potential. 
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5 Fermionic Determinant w.r.t. Boundary Conditions 



In this section we shall compute the ^-dependence of the fermionic determinant. We shall 
see that the scattering of the fermions off the boundary generates a CP-odd 0-term in the 
effective action for the gauge-bosons. 

5.1 Zetafunction Definition 

The Dirac operator and the boundary conditions are both flavour neutral. Thus the deter- 
minant is the same for all flavours and it is sufficient to calculate it for one flavour. For the 
explicit calculations we shall use the gauge-invariant (^-function definition of the determinant 



mm 

logdete(i^) := ^logdete(-^) := -~ 
2 2 as 



Ce(-^,s) (32) 

s=0 

and calculate the ^-dependence of the C-function by means of the boundary Feynman Hell- 
mann formula (|9|). Denoting {/Xfc|/c € N} the (positive) eigenvalues of — the corresponding 
C-function is defined and rewritten as a Mellin transform in the usual way 

oo 

Ce{s) :^ G(-^,.s) :^ E/^r = f7^/ 1^"' tr ,(e-*(-^)) dt (33) 
for Re(s) > d/2 = 1 and its analytic continuation to Re(s) < 1. 



5.2 Stepwise Calculation 

We will study how det0{ilpa,c) varies with 6, a and c to compute the normalized determinant 

dete(i^) _ dete{i]pa=i,c) 



deto(«^) deto(i^a=o,o) 
The calculation is done in three steps. We shall calculate all three factors in 

dete(i^ct=i,c) ^ dete{ipa=i,c) _ deto(z-^a=i,c) _ deto(i-^a=o,c) 
deto(i-^Q:=o,o) deto{ipa=i,c) deto(i-^a=o,c) deto(i-^Q=o,o) 



(34) 



(35) 



in turn. 

From the generalized Feynman-Hellmann formula (|9|) and the fact that ip has no zeromodes 
so that the various partial integrations are justified, the ^-variation of (|3^) is found to be 



oo 



^Ce{s) = ^ f-' tre(e*^75). (36) 







de^"' ' r(s) 

Now one can use the asymptotic small-t-expansion for e*^ /, where / is a testfunction, 

tre(e*^/) = ^ E i'"^' tr//a„/2(/)+ /'6„/2(/)) , (37) 

m=0,l, 



and where the am/2i^m/2 denote the corresponding volume and boundary Seeley DeWitt 
coefficients respectively. Plugging this into the expression (^) yields |jl^, 21] 

^ilogdet,(-^) = -i-| tr(ai(75))-^/tr (61(75)). (38) 
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For the squared Dirac operator —Ij^ that part of oi which leads to a nonvanishing 75-trace is 
known [22| to be eFoi/27r, i.e. independent of 9. On the other hand §hi{.), which depends 
on the boundary conditions, is calculated explicitly in the appendix to be 



2 i y 1 y sinh(6l) y -1 e 

and does not contribute for / = 75 
Integrating with respect to 9 yields the following first factor in (|3 



'a=l,cJ 



exp{-|^ / = "^""P^"^ / 



deto(i^a=i,c) 

To find the a-variation leading to the second factor we use dl^) in (|33|) with the result 



00 



^G(-^«,c,«) = y f-4r,(e*^"2<A75), (41) 

where we integrated by parts. Again we use the small-t-expansion (37) of the heat kernel, 
but now with test function / = 2^75. Thus 

^llogdet,(-^„,J = tr (ai(2075)) + i- ^ tr (61(2075)) (42) 

where the universal ai(.) yields the wellknown Schwinger term [p^]. Since now the normal 
derivative of the testfunction on the boundary in non-zero, the last surface term contributes. 
Using ( |39[ ) we end up with 

^ilogdet,(-^„,J = i-| 2c/>A<^- i- j log(e^) a„</>. (43) 

M dM 

Setting = and integrating with respect to a yields the following second factor in (35): 

deto(.^.=i,e) ^^_l;^a^_ (44) 
deto(?-^Q=o,c) 

We are left with the task to calculate the third factor 



deto(i-R,=o,c) I f d a ^ 
deto(^^a=o,o) 2 J ds s=odc 

For that we computed the heatkernel of the operator 

-PLo,c = -{{do - 27ric/(3f + d!)l2 
for 9 = 0. The explicit result is 

K{t,x,y) = ^ ^ (_i)™+"e-(?°-"/3)V4V.ic(cO-„/3)//3 



zxz 



g-(^i-2mL)2/4t _g-(»;-2mL)2/4t 
_ -{»?+2mL)2/4t „-{Ci+2mL)V4t 



(46) 
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which results in the trace (V = I3L) 



/V / 1 — , (n/3)2 + (2mL)2 

tr,=o {K{t,x,x)) = ^(l + ^(-ir+"e- ' 4^ cos(27rnc: 

M 



(47) 



where the prime denotes the omission of the (m, n) = (0, 0) term. The c-derivative of the 
Melhn transform, after substituting t 1/t, reads 



V 



27rr(s) 

u 

which may be integrated by parts (for s > 0) to give 



J ^'(_i)-+ng-t[{n/3/2)2 + (™L)2] d_ ^^^^^^^^^ ^-s (^g) 



d_ 

dc 



2V s 
'Vf(s) 



"^+"^ — cos(2^nc) t-'-^ dt . 



[(n/?)2 + 4(mL)2] dc 



Only the pole of order one of the integral at s = can contribute to the s-derivative at s = 
of the ("-function. Since this pole entirely stems from the lower limit of the integral we may 
split the latter into two parts 



d 



dc 



to obtain 



rCo(-^=0,c>^iO) 



2V s 



TT r(s) Jo 



... + 



2V 



— (s+w +...)■ 



m+n dc 



cos(27rnc) 



[(n/3)2 + A{mLf 



d 

ds 



Plugging this result into (^) we end up with the expression 



- r(c) = log 



deto(i^a=o,c 
deto(i-^c 



'a=0,Oj 



Y_ v7 , ym+n cos(27rnc) - 1 
' TT ^ ^ ^ (n/3)2 + (2mL)2 



(49) 



(50) 



for the third factor of the functional determinant (|34| ) in the factorization (|35D . With the 
help of 



Im(r) / e2'^*("*'^i+"''2) 



vr 



m + rn ■ 



-2 log 



1 



ry(r) 



2+01 

^2 



this result can be rewritten as |2S] 

-r(c) _ deto(i^a=o^c) 



deto(i^Q,=o,o) 



9z{c,iT) 
ez{0,ir) 



(51) 



These two equivalent forms will be useful in the low- and high- temperature expansion of the 
condensates. 
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5.3 Effective Action 

Now we can combine the classical (euclidean) action of the photon field, rewritten in the new 
variables ( [T^ ) 

\f^uF^, = ^A0A(/> = : Sb^ (52) 
with our explicit result for the functional determinant p^). Collecting the contributions 



(|40|,^,^l|) and adding the classical action (^) we end up with the effective action 

r = re[c,0] =iV/r(c) + reM (53) 

where r(c) has been given in ( |5l| ) and Tg[(p] is 

Te[4>]^^y ^A'4>-fi^J 4>^<P + e-fi^ J A(t>^ (54) 



MM M 



and 



:^ (55) 

is the analog of the ry'-mass in QCD. We have used that the functional determinant is the 
same for all flavours. The functional measure takes the form 

dfie[A] = ^ e-r«['=.'^] dc D<p 5{x) Dx • (56) 

We have taken into account that the gauge-variation of the Lorentz gauge-condition 

F ■= d^A^' = Ax 

and the Jacobian of the transformation from {A} to the new variables {(j), c, x} are indepen- 
dent of the fields. Actually, the corresponding determinants cancel each other. 

We conclude that the expectation value of any gauge-invariant operator O (which will not 
depend on x) is given by 

^ '~ JdcDcp e-r«M] ^^^^ 



with rg[c, (/)] from (53,54). 



6 Chiral Condensates 

Our result ( |57|) may be applied to calculate the chiral condensates as 

(^t(,)p^^(,)) = j^,^^ ,-r,[,,] (58) 



with Se from ( |3lD and from (pg). Both the (exponentiated) action and the Green's function 
factorize into parts which only depend on c and on (p, respectively. Thus ( p8| ) factorizes as 

{ipHx)P±ij{x)) = C^{x) ■ D^{x) (59) 

with x'^-independent factors 

, , _ JdcSe{x,x)^ e-^/^(-) , , _ / e^2^(-)~r.M 

^ (^)- e-^/rW ' ^ ) - / Dcp e'^^M ^^^^ 

which depend on the parameters 9, Nf, P, L. Here and below the c-integrals extend over one 
period, e.g. [-1/2,1/2]. 
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6.1 Harmonic Integral 

Now we shall see, how far we can evaluate the first factor in (|5^) which contains the integrals 
over the harmonic part of the gauge-field. 

Plugging in the Green's function (p5| , ^ ) as well as (|5^) we obtain the unevaluated expres- 
sion 



1/2 „, ,fc + ; coB(27ric)-l 

/ cos(27rnc)e ^'^ ^ ^ ' dc 



AttL (,^ + m)2 + (nT)2 1/2 y^', sfc^; cos(2,ric)-i 

-1/2 

(61) 

To investigate the low-temperature expansion we use ( p9[ ) and the upper line in ( |5l] ) and 
arrive at 

4L sin(vr[e - ^nr]) ' j </(c,zr) ' ^ ^ 



Alternatively, for the high-temperature expansion we use (|30| ) and the lower line in (plD, so 
that 



2/5 sinh(^[C + m]/T) / dce-^'^//^ <^(ic/r, i/r) 



(-1)™ J g-7rc[7V^cT2(C+m)]/r ^^^/ (^g/r, i/r) 



For one flavour the c-integral in is easily calculated and one finds 



P^S ( _-\\n p-TTT'n? 

cHx^) = ^-TrY. S i ■ n • (64) 

sm vr f — inr 



6.2 Nonharmonic Integral 



Now we shall compute the second factor in (59) as defined in (|60|). We recall that the 
integration extends over fields 0, which are periodic in the and satisfy Dirichlet boundary 
conditions at the ends of the cylinder, i.e. atx^ = 0,L. 
Doing the gaussian integrals one ends up with 

27r 

~f 

where the integration is over x' and the kernel 



D^{x^) = exp{^i^^2(x,x)}-exp{±^.y' A' K^2{x,x') ± ^ ■ J A'K^2{x',x)} (65) 



^2 11 

K,,2(x,y) = (x\ — —. — ktIu) = (x\ — —\y) — (x\ — kIu) (66) 

is with respect to Dirichlet boundary conditions. Being the difference of two Green's functions 
with the same singular behaviour it is finite at coinciding arguments. 
The explicit form of the kernel is 

K^,2{x,y) = :^J2m,n&z ((2mL)2 + (n/?)2 " {2mLy + {nf3y + {fiV/Try) ' 



cos 



(— ^ j (-^ j sm , (67) 
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where the prime indicates the omission of the term with m = n = 0. For coinciding arguments 
K becomes x'^-independent as required by translational invariance. For performing either the 
sum over m or over n in (67) one uses the formula 

cos(jx) TT cosh(a(7r - x)) ^ rn o-^n 

+ a sinh(a7r) 

to end up either with the expression 

(x, ^) = ^ E (^^^^^ - (n ^ ^n^ + (/xL/7r)2)) (l - cos(27rnO) , (68) 

which is useful for the low temperature expansion, or alternatively with the expression 

/ N 1 coshfmTr/r) - cosh(m7r(l - 20/t) , / ' , „ ,„ 

^ 27r m smh mTT/r 

m>l ^ ' ' 

, cosh(ML(l-2g))-cosh(/.L) 
2r ^ 2;U/3 sinh(^L) ' ^ ' 



which is useful for the high temperature expansion. Both expressions (^8|) and ( |69| ) do indeed 
vanish as x^ reaches the boundary in accordance with the imposed boundary conditions. 

Once we have the explicit formula (|67| ) at hand we can compute in a straightforward way 
the expression 



/4 . 1 n 

A,A>(...)<i. = -- j: (---^-^-^)si„(™0, (70) 



4 ^ ,1 

?i=l,3,.. 

Applying the formula 



En sin(nx) vr shfafvr — x)) + sh(ax) , -, r n n 

2^2 = TT^ X G 0, ^ 71 

n=l,3,... ^ ' 

the expression ([70|) is seen to take the simple form 

A jv- . N , sinh(/xL(l - ^)) + sinh(^Lg) 

A^iC 2(z,x) = . , . ^ • v^^) 

^ smh(^Lj 

6.3 Final Result 

Now all pieces to compute the chiral condensate (|5^ ) have been calculated. For we have 
the two alternative forms (62) and (^), and is given by ( |65| ) wherein we can use one of 
the equivalent representations (|6|) or (^9|) for K^2 together with (l?^) . Thus we have 



" ^ =^^(^1^ - ^^^]) / e^^' (c, 



-pf ]^ E (^^!^ - (« - + (^)=)) (i - ■=«(2-»{)) ) 
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smh(^L(l - g)) + sinh(//LO 

exp{±6' ^TT-Ta J v73) 

smn(/iL) 



1 (-1)"^ /dce-^^[^/^^2{C+m)]/r ^^^/(^^/^^^/^) 



, 1 coshfrnvr/r) - coshfmvrfl — 20/t) , / „ //^/Snom 

expiiT^ / ^-n r\ ("^ ^ \/"^ +(7;— ) )} • 

Nf m sinh(m7r/r) ^ V 27r^ 

^27r cosh(^L(l-2g))-cosh(A.L) x^ 
""""P^iVjV 2r ^ 2;u/3 smh(/iL) 

smn(/iL) 

with excellent convergence properties in the low- and high-temperature regime, respectively. 

This result is one of the two main results of this article. To simplify the analysis we shall 
now study the condensates at the midpoints of the cylinder. 

6.4 {ij^P±ij) at Midpoints 

If a condensate survives at the midpoints when the boundaries are taken to infinity then the 
chiral symmetry is broken. 

For = L/2 the formulas (|6^), (|6^ ) simplify to 



r±(^\ - I1 I oY^r 1^n /cos(27rnc) 9^f{c,iT)dc 

^ \ cosh(n7rTj J 0^ [c,iT) dc 

2 l^nNf 



(75) 



C^{k) = ±^^V( i)»^ /^osh((2m+l)7rc/r)e-^/^/-g3''^(ic/r,z/r)dc 
2 /3 ^0 sinh((2m+l)7r/2r)/e-"^/^'/^0f^(ic/r,Vr)dc' 



The formulas (pq) and (p9) simplify to 



■2' 27r V sh(m7r/T) ' V '27r' 7 8r 2///3 sh(^L) 

Depending on which one of the equivalent forms ( |68[ ) and ( |69| ) for K^2 on the diagonal is used 
the factor at the midpoints is found to read 

2 '-A'f "7-; n V 71" J 

J n=l,3,... ' 

exp{ - l/ch(^L/2))} (79) 

n±/^^ f 1 ^ th(m7r/2T) /TTT'^^ \\ 



which can be used to derive the low and high temperature expansions, respectively. 
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7 Noncommutativity of the Limits (3 ^ ^ and L ^ oo 

In this section we show that the condensates at the midpoints, P±il>)i3{^), depend on the 
order in which the hmits /3 — > oo and L ^ oo are taken; for Nj = 1 the condensates survive 
only if we first let (3 —> oo. 

7.1 Limit (3 ^ oo for finite spatial length L 

Here we derive the low temperature limit, i.e. the condensates for /3 large compared to the 
fixed spatial length L and the induced mass /i. 

From the explicit expression (^) we see at once that 

T 

C^q) = ^-^i^ + 0{e-^-''^'^'^)) (81) 

for any number of flavours. 

In order to get the corresponding limit for the second factor D^{^) in (^) we use ( ^9|) 
and perform the asymptotic expansion of the coth to get 



exp{ =f6'(1 - l/ch(/iL/2))} (82) 



adapted to /3 3> L as an intermediate result. Using the identity |27] 



valid for x > 0, where 7 denotes the Euler Masceroni constant and Kq the zeroth Bessel 
function the second factor can be rewritten as 



D^i^) = e^/^/(^L/vr)'^'^%xp{- J-^(-l)^Ko(i/xi.)} 



exp{T0(l-l/ch(^L/2))}.O(exp(— e-2-/^/2^)) . (84) 
Combining (^) and ( |84|) we get the result 

{i^'P±i^),i^) = ±^(f^)^/^^e^/^/exp{-|-^(-l)^i^o(j^L)} . 

exp { ± ^/ch(^L/2)} • (1 + 0(e-2-/3/2i)) (35) 

where the 9 dependencies are found to cancel up to exponentially small remainders. In 
particular we have found a nonzero value for P±ip) for midpoints at zero temperature for 
any Nf for finite spatial length L. 
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7.2 Limit L ^ oo for finite temperature 

Here we give the large volume expansion of (^9|) valid for length L which is large as compared 
to the fixed inverse temperature /3 and /j,^^. 

The first task is to derive the high temperature asymptotics for the first factor in (|59|). 
By a variety of manipulations including infinite product representations for the exponential 
factors constituting the measure we arrived at the asymptotic result [24] 

{Nf = 1) 



C^(-) = <j ±'—-2e 



±V.4e 



2]Vf-l „r 
-2 i — 



{Nf = 2) 



(86) 



/ " {Nf > 3) 



which is an exponential decay which goes faster as the number of flavours increases. 

Also, we performed the asymptotic expansions of the hyperbolic functions in (^) and 
arrived at 



/ 



m 



/ m>l 1>1 

1 + 2 E(-l)'e-'^^ 
r 1 vr / i>i 



/ rmr/T -«7r^m2+(/i/3/2)2/7 



/xL/2 



)}■ 



exp{:f 0±20^(-l)'e- 



(2^+l)^tL/2^ 



(87) 



where everything is at least exponentially suppressed as compared to the growing factor in 
the second-last line. 

Combining (^) and (^^ we end up with the result 

5 TvL 

2— {Nf = 1) 



1 



{Nf = 2) 



Try/2(3L 

l-e~'^~ {Nf>3) J 



-p{-^E^o(././?)}-0(exp{-A,-W/3} 



exp • 



/^-(2«+1)mL/2| 



which gives a decay 



ibconst • -^e (A^/ = 1 



± const • e 



{Nf > 2) 
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for spatial lengths L which are large compared to the inverse temperature (3 and the inverse 
charge e~^. 

7.3 Noncommutativity of the limits (3 ^ oo and L ^ oo 

Using the results of the previous subsections it is easy to show that the two limits /S ^ oo 
and L ^ oo do not commute. 

Recall that {ip^ P^ip) {^) is a shorthand for {i^^ P±il^)e,Nf,(3,L{x^ '■= ■§)• 
Now the formulas (||), (||) imply 



±^e^JJ^ {Nf = l) 
hm hm (^tp^^)( ) = <^ " " (90) 

lo {Nf>2) 
lim lim ('0"^P±V)(-) = ( V iV/ > 1) (91) 

respectively, which is the other main result of this paper. 

From a physical point of view this means that the system under consideration shows a 
distinctive hysteresis phenomenon: When both of /? and L are sent to infinity, the one-flavour 
system keeps the knowledge of which limit was performed first in the actual value of its chiral 
condensate. Obviously there is no such non-commutativity for finite changes of the lengths 
P and L. We shall further comment on this interesting behaviour in the conclusions. 



8 Discussion and Conclusions 

In this paper we have performed in a functional framework the quantization of the Nf flavour 
euclidean Schwinger model inside a finite temperature cylinder with SU{Nf)A breaking local 
boundary conditions at the two spatial ends to trigger chiral symmetry breaking. We have 
determined the effective action for the bosonic subsystem subject to these boundary condi- 
tions, which arises after integrating out the fermions. We have shown the way the expectation 
value of an arbitrary gauge-invariant operator can be computed and in particular we have 
performed the calculation of the condensates P±ip){x) (to be used as the most simple order 
parameters) for any point x inside the cylinder and any value of the inverse temperature /3 
and spatial length L. 

The quantization was greatly simplified by the fact that the boundary conditions chosen 
(the CB-boundary-conditions) completely ban the zero modes. Once more we emphasize the 
fact that our results have been obtained purely analytically and without doing 'instanton 
physics'. The technical aspects are rather different as those one encounters when quantizing 
the theory on a sphere |l^ or on a torus |lj, . 

Nevertheless our results are in full agreement with the earlier instanton-type and small- 
quark-mass calculations. Thus it seems that the CB-boundary-conditions applied at the two 
spatial ends of the cylinder give a perfect substitute for introducing small quark masses to 
trigger the chiral symmetry breaking and a real alternative to the study of torons [^] or 
fractons P or singular gauge- fields on [13|. The real advantage is of course the fact 
that they constitute almost exactly the border of what can be calculated analytically. The 
functional integral over the prepotential is gaussian, whereas, in general, the integration over 
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the harmonic part of the gauge-potential is not. However the latter reduces to gaussian 
integrals in the low and high temperature expansions. 

In the low temperature limit = l/y^iVje^/vr <^ L <^ = T^^ we found for the chiral 
condensate the asymptotic value 



which, when restricted to the two-flavour case reduces to 



2' V 16vrL 

This expression is identical to the result of Shifman and Smilga [P, who allowed for fracton 
configurations on the torus. 



In the high temperature limit T = (3 ^ ^ ^Njc'^/tt ^ L ^ we found for the chiral 
condensate an exponential decay with T. 

For intermediate temperatures T = (3^^ ^ J Nje^ /t: and finite L one has to retreat 



to numerical methods to evaluate the remaining sum and the integrals in (pBI) and (79) or 



equivalently in ( |76[) and (^). One realizes that the observable {4>P±il^) viewed as a function of 
T strongly resembles the behaviour of an order parameter in a system which suffers a second 
order phase transition for the case Nf > 2. However, the chiral condensate does not really 
vanish at any finite temperature, it is just exponentially close to zero for temperatures larger 



than the induced mass = ^NfC^ /tt. Thus, in a strict sense, the chiral symmetry remains 
broken even for Nf > 2 at all finite temperatures as long as L stays finite, as has been argued 
to be a general fact by Dolan and Jackiw Q. However, if L is sent to infinity for finite /?, the 
condensate exponentially drops to zero. 

Our main result is the fact that the limits /? — > oo and L ^ oo do not commute for the 
observable {ip'^P±ip) in the Nf = l case, since 

hm hm (V^tp^^)(^) = ±^e^\ - {Nf = l) (94) 

L— >oo /3^oo Z 47r y TT 

lim lim (V^1'p^^)(:^) = q ( V iV/ > 1) (95) 

/3— >oo L— >oo 2 

which implies that there is no unique infinite volume limit. Thus it seems that the combination 
of finite-temperature and CB- boundary conditions provides an interesting tool for driving 
this system either into the true or the wrong vacuum state. The result (94,9^) is rather 



remarkable, since it means that the one-flavour system shows some hysteresis phenomenon: 
As far as we are aware of the literature, such phenomena are known for spin systems but they 
are rather untypical for analytically solvable field theories. However one of the interesting 
new results in this respect is the work by Hetrick, Hosotani and Iso about the massive multi- 
fiavour Schwinger model on the zero temperature cylinder [^]. They analyzed the situation 
for small quark masses and finite (cyclic) spatial length L. In particular they found that the 
two limits m — > and L — > oo fail to commute. Thus we conclude that chirality breaking 
boundary conditions give an interesting alternative to introducing small quark masses. 
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A Explicit Construction of the Fermionic Heat Kernel 

In this appendix we sketch the construction of the heat kernel of the squared Dirac operator 
(^-Pla=o)^ = + 27rc//57^)^ on a thermal manifold, which allows to compute the relevant 
Seeley DeWitt coefficient, a task, which itself is postponed to appendix B. For that we con- 
struct the heat kernel K on the finite cylinder {{x^,x^) | G [0, /?[ , 2;^ > } which obeys 
(besides the usual heat kernel relations) the boundary conditions 

(5,K)(t,x°,0,y) = i^(t,x°,0,y) (96) 
{Beim{t,x\Q,y) = (i^K)(t,a;0,0,y) (97) 
k{t,x^+P,x^,y) = -K{t,x^,x\y) (98) 

as well as the adjoint relations with respect to y, where Bg is a shorthand for Bi{6) defined 
in (|10|). 

The trick is to start considerations on the half plane {{x^ , x^) \ > 0}, since here the above 
squared Dirac operator can be decomposed as 

(i^^ + 27rc//3 • a^f = e^^'"'"' {i^^fe-^^'"'"' (99) 

and correspondingly the free heat kernel takes the simple form 

_}_f.-{{e? + {e?)/it^2^ice/l3 ^ J_g-{(CO-47r»ct//3)2+(5i)2)/4tg-4^2^2i//3' (100) 

where = x^ — y^ ,£} = x^ —y^. Using that the kernel can be Fourier transformed and from 
the mirror principle one is led to consider the expression 



1 



(2vr) 



2 



oo oo / , 



^(2^)2 J J \ giko,ki) h{ko,ki) 

— oo— oo ^ 

as an ansatz for the heat kernel of the operator (i^)^ = —A • I2 on the half plane. The 
boundary condition at = immediately transforms into an algebraic relation among /, g, h 
which is solved by the expressions 

9{ko,ki] 
h{ko,ki] 



e'^^{kQ-iki) 


- (ko-iki) 


e^S{ko+iki) 


- (ko-iki) 




iki 


e-^\ko+iki) 


- (ko-iki) 


e^\kQ+iki) 


- (ko+iki) 



e'^^{ko+iki) - {ko-iki] 
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The resulting integrals can be done in two steps. First only the numerators of the functions 
/, 5, h are taken into account and the resulting expressions are integrated over. Second the 
full expressions have to be read as differential equations in x^,x^ in the manner indicated by 
the previously omitted denominators of the functions /, g, h. There is a unique solution to 
this procedure which falls off in both plus the positive x^ directions (note G R) : 

e-^^^^^^%^.(l + erf(^^^%^ 

Since on the half-plane the operator (i^ + 27rccr^//3)^ has the decomposition ( p9| ) this im- 
mediately yields its heat kernel ( to be denoted K ) hy just including a factor e^'^*'^ in 
each term. Finally the finite temperature boundary condition ( p8[ ) is taken into account by 
substituting by SP—n(5, including an additional (—1)" and performing the sum over n ^7^. 

The heat kernel Ka=Q of {ip\a=oY — (^^ + 27rc(T^//3)^ subject to the boundary conditions 
(PI) " (PI) on the half cylinder { (x°,x^) | x'^ G [0, /?[ , > } takes the final form 



e 4t g27ric(eO-n/3)//3 



4:TTt I -ch6' -e-'^sh6' 

((e°-nfflch^+,^ sh0f / ^J(^O-n/3)sh0-r?ch0A , , 

e • ( 1 + erf(^^ "^tU^ V 

where the sums run over n£Z and can be seen to converge absolutely and thus uniformely. 



B Extraction of the Relevant Heat Kernel Coefficients 

In this appendix we shall compute the surface Seeley DeWitt coefficient bi of the operator 
— which enters the calculation of it's functional determinant. We first note that in general 
§ tr (bm{f)) with a smooth test function (p on a d dimensional manifold M has the expansion 

j tr {bM) = T. {bm.p{R, X, F.) ■ dl^) , 
dM P=^dM 

where b^.p is a gauge-invariant and Lorentz-covariant local polynomial in the intrinsic and 
extrinsic curvatures of the boundary as well as in the field strength and its covariant derivatives 
on the boundary. Here d^ip denotes the p fold derivative of the test function 99 along the 



19 



(outward oriented) normal of the boundary. In the case of a two dimensional manifold with 
Hrasko Balog boundary conditions ([lO|) the expansion of / tr {bi{ip)) simplifies to 



tr (6i(^)) = ^ tr (6i.o(^) x ■ + f tr • dnf) . 

For our purposes it is sufficient to know the coefficient bi i, since the first term does not 
contribute to (42) (due to (p := H+H^ = on dM) and in (^) it would yield an uninteresting 
constant which finally cancels in expectation values of gauge-invariant operators. The function 
61,1 can be determined from the heat kernel on the diagonal, K{t,x,x) of —fP = —I/9'\a=i 
which is identical to K{t,x^x) of ~IJP'\a=o by calculating 

00 

j K{t, X, x) ■ ip{x) = j k{t, X, x) ■ ip{x) ~ J kit, X, x) ■ (((3(x°, 0) + x^ ■ diip{x°, 0) + dx^ 

M M 



where K denotes the heat kernel ( 101 ) calculated in appendix A. In writing this expansion we 
have anticipated that for small t the heat kernel on the diagonal is sharply peaked about the 
boundary whereupon it is justified to expand the test function (/? about x^ = 0. Using this 
result and denoting ip'{x^, .) the first derivative of 93 with respect to it's second argument one 
has to compute an expression whose first few terms in the small t expansion take the form 

/ 



e^sh — ch \ 1 





00 



, -ch -e "^sh / A-Kt 





00 

.2 , 



^ ( e^'sh -sh \ /■ chn/3 - 2«shx (chnl3~2is\YX)^ 

+ Z^(-1) 1 _sh +e-^sh j ■ J 87rV2if3/2 ^ 



/ ch2x + ishn/3 \ _2mnc 7 / 



+ E(-i) 



e^sh -sh \ [ chn/3 - 2ishx (ch?T,/3-2zsha;) 



-sh +e-'^sh / J 87ri/2if3/2 ^ 



2 



(l-erf( ^^'^^+;f^^ )) dx.,.'(x°,0) 



where the first line gives the usual oq coefficient whereas the remaining four integrals contain 
information about the 61/2 and 61 coefficients. Here and below we use the abbreviations 
sh = sh6', ch = ch6'. 

The first and second integrals are easily evaluated using the formulas 

00 , 00 



Ii := / e t dx = — , I2 '■= I X ■ e t dx ^ 
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The third and fourth integrals are handled using the formulas 
h: = 



chn/3-2ishx (chn/3-2zsha:)^ ^.ch2x + ishnjS ^. , 
e 4t (l-erf( )) dx 



87ri/2it3/2 



1 ch "^/3^ 1 1 _chf_ji^ fshn/3 

e 4' + ,^ ,^ — • e 4t eric(- 



87rV2ti/2 sh 



87ri/2ti/2 sh 



2ti/2 



chn/3 - 2ishx _ {chnf]~2ishx)^ ch2x + ishn/3 , , , 
X ^ wo. o/o e 4t (l-erf( — ^7;^ 



87ri/2it3/2 



1 ch 1 1 

p 4t -I 

87r sh 87rV2ii/2 gh 



2*1/2 

(chnP~2ishx)^ ch2x + ishn/3 



"erfc(- 



2ii/2 



) dx 



which result in the small t asymptotics 



I3 = < 



1 ch -ILiZ_ i 

87ri/2ji/2 sh ^ 47rsh' 

1 ch I 1 1 

'8^i/2jl/2 sh 8771/2*1/2 sh 



^e- — (l + 0(t)) (n>0) 

(n = 0) 



' ^'^e-'^ + i^e-4f^(l + 0(t)) (n<0) 



'87ri/2tl/2 sh 



(102) 



h = { 



87rsh ^ 



87rV2sh2n/3 



e it 



■(1 + 0(^1/2)) (n>0) 



1 ch , 1 log(ch+sh)~log(ch-sh) 



" 87r sh Btt 



2sh^ 



STTsh ^ 87rl/2sh2„/3 



(n = 0) (103) 



J^e-4f^(l + 0(tV2)) (,<o) 



where the result for immediately follows from the asymptotic expansion pSl 



z e^'erfc(z) - 1 + f^i-l)''^-^--^^^^! (^ ^ 00, |argz| < ^) (104) 
fe=i ^ ' 

whereas the expression for /4 results from a computation establishing the asymptotic be- 
haviour 



f{w) = y e-(=^'^-^^^^')'erfc(chx + ish'u;)dx 


_ -w^ f i I ch 1 ill 



(105) 



for i« » 1. 

Putting everything together we arrive at the small t expansion of the heat kernel 



K{t,x,x)ip{x)dx -Oit^/"^) 



(106) 



M 
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+ 4^ • I E(-l)''e — — cos(27rnc)| • / x^) d^x 

+ 5^I7W • { ( V -1 ) cos(27rnc) + ( ^"i ) } • / dx' 

'i -e E (-l)"e"— sm(27rnc) ^ -/^(x^O) dx" 



which is used to determine the effective action (|53D. This formula resolves also the appar- 
ent paradox that the 0-term in the effective action ( |53D is linear, whereas the whole model 
was defined through hyperbolic functions of ^, thus there must be an invariance under the 
replacement 9^9 + Ini. 
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